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A molecular dynamics simulation is performed to investigate spatial scale of low energy excitation 
(LEE) in a single linear chain of united atoms. The self part of the dynamic structure function, 
Ss(<7,w), is obtained in a wide range in frequency space (uj) and reciprocal space (q). A broad peak 
corresponding to the LEE is detected at u/2-k — 2.5 x 10 11 s (= ojlee/27t) on the contour maps of 
Ss(q,u>), near and below the glass transition temperature (T g =230 K). The Ss(g, wlbe) is symmetric 
around a maximum along the logarithm of q. The inverse of g ma x, giving the maximum position of 



Ss(5, wlee), depends on temperature as 2n/q n 



T u " 2 for 60 K < T < T g and 2n/q n 



for T g < T < 600 K, which is the spatial scale of the motion corresponding to the LEE at low 
temperatures. Based on a Gaussian approximation for the displacements of monomer groups which 
give rise to the motion relevant to the LEE, it is found that the number of monomers contained in 
a group is about 6. 

PACS numbers: 61.20.Lc, 61.41.+e, 61.43.Fs, 83.10.Rs 
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Low energy excitation (LEE) is observed universally in 
various glassy materials, such as amorphous metals 



network glasses 
glassy polymers 




molecular glasses (ESQ' an d 
jjj. It has been studied inten- 



sively with the relation to anomalous thermal properties 
of the glassy materials in the temperature range from 10 
to 30 K. 

In many glassy materials, we can observe LEE as a 
broad peak in neutron and Raman scattering data at low 
energy range from 1 to 5 meV near and below the glass 
transition temperature (Tg). Because the temperature 
dependence of the peak height can be well fitted by Bose 
factor, it is often called as "Boson peak". LEE can be 
recognized as an excess excitation superposed on that due 
to the vibrational density of state for a crystal (Debye 
solid), which is observed as an excess specific heat in the 
temperature range from 10 to 30 K LEE affects the 
thermal conductivity in the temperature range, because 
the coupling between phonons carrying the heat and LEE 
occurs when their characteristic frequencies are close [l2l . 

IB 

So far, several models have been proposed for the ele- 
mentary processes of LEE in order to explain the unusual 
behavior of the specific heat and thermal conductivity in 
the temperature range from 10 to 30 K 0, M M- Ex- 
perimental researches were performed to obtain the char- 
acteristic sizes of LEE for various materials 0, E3j Ell 
using the proposed models. Although many results have 
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been reported, the evidences justifying the models are 
not sufficient. 

Recently, much attention has been paid to the origin of 
LEE, because significant improvements of experimental 
techniques enable us to examine the scattering intensity 
in a wide range of energy- momentum space [191 ] . It im- 
plies that we can in principle identify the characteristics 
of LEE directly without any specific models, although 
the unified understanding of the origin of LEE has not 
yet been attained. 

On the other hand, the energy and momentum ranges 
of LEE are accessible by molecular dynamics (MD) sim- 
ulation. Generally, MD simulation is a powerful tool 
to investigate atomic scale non-equilibrium phenomena. 
Thus, it is expected that the problems of LEE can suc- 
cessfully be clarified by MD simulation. Stimulated by 
the models proposed and the experimental results, MD 
simulations for various systems have been performed 
(20l l2ll [22I [23l | , where suggestive results have been re- 
ported. 

At the present stage, the information on the character- 
istics of LEE independent of the models is highly required 
in order to elucidate the nature of LEE. The purpose of 
this work is to examine the spatial scale of LEE by MD 
simulation. We calculate self part of dynamic structure 
function which can be compared with experiments di- 
rectly. We choose a polymer model system to investigate 
LEE, because polymer is one of the common materials 
which easily vitrify, and its LEE has been studied exten- 
sively. 

We employ a united atom polyethylene (PE) model 
as the molecular model for the simulations, which is 
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essentially the same as that of our previous work [24|. 
We consider energies of bond stretching, angle bending, 
and torsion of the model PE. Truncated Lennard-Jones 
12-6 potential is used to reproduce van der Waals in- 
teractions between two united atoms apart more than 
three bonds, where a cut off radius is set to 1 nm. The 
bulk polymer system is made of a single linear PE com- 
posed of 2000 united atoms placed in an MD cell under 
three-dimensional periodic boundary conditions. New- 
ton's equations of motion are integrated by leap-flog al- 
gorithm with a time step of 4 fs which is 1/20 of a period 
of the bond oscillation. The temperature and the pres- 
sure are controlled every 1 ps by ad hoc velocity rescal- 
ing method and loose-coupling method [25[ , respectively, 
where the correction for the pressure due to the cut off 
is taken into account. 

An initial configuration of the model chain, with fixed 
bond lengths and bond angles, is generated using tech- 
nique developed by Theodrou and Suter 26], and then 
a short MD simulation of 10 ps at a constant volume 
(1.45 cm 3 g _1 ) and 600 K is performed to reduce the 
overlaps of the atoms. This state is then fully relaxed for 

1 ns at 1 atm and 600 K. The final state of 600 K obtained 
above is quenched to various temperatures between 200 
K to 500 K, and the simulations of 2 ns are performed at 
a constant pressure (1 atm) and the constant tempera- 
tures. For the simulations at the temperatures below 200 
K, the final state of 600 K is quenched to 200 K, relaxed 
for 1 ns, and then quenched to the temperatures between 
60 K to 180 K again. After the stepwise quenching, the 
simulations of 2 ns are carried out at a constant pressure 
(1 atm) and the temperatures. We adopt the stepwise 
quenching procedure at the low temperatures to repro- 
duce the temperature dependence of the specific volume 
of amorphous PE [24[ . 

Here we are interested in the LEE in the isotropic and 
homogeneous amorphous state of PE before it crystal- 
lizes, so that the total time of the simulation must be 
shorter than the characteristic time of the crystalliza- 
tion. In the present study, the total simulation time of 

2 ns at each temperature is much shorter than the time 
scale of the crystallization which is about the order of 10 
ns observed previously [U HI] . 

In order to investigate the spatial scale of LEE, we 
evaluate self part of dynamic structure function, which 
can be measured by various scattering experiments. The 
self part of the dynamic structure function is defined as 



/+oo 
F s (q,t) e-^dt 
-oo 



(1) 



where q is the wave vector, to is the angular frequency, t 
is the time, and i*s(q, t) is the self part of intermediate 
scattering function defined by 

1 N 

F s (q,t) = — < £V iq(r ' (t+to) - r ' (to)) > t0 (2) 
where ./V is the number of the united atoms, Tj(t) is the 




FIG. 1: Ss(q, to) vs. lo/2-k. The results of the temperatures 
from 200 to 460 K are plotted. The arrow indicates the posi- 
tion of the LEE. 



position vector of jth united atom, and < ... > to denotes 
the average about t . In our study, we employ an algo- 
rithm developed by Matsui et al. [29| to obtain 5s (q, to), 
where the 5s (q, to) is calculated during the final 1 ns at 
a given temperature. After calculating the 5s (q, u>), we 
take an average of it about the solid angle of q, and 
then obtain Ss(q,u>), where q = |q|. We set the q and to 
ranges to be 5.3 x 10 s m _1 < g/27r < 4.6 x 10 10 m _1 and 
9.5 x 10 8 s- 1 < uj/2tt < 2.5 x 10 14 s" 1 (3.9 x 10" 3 meV < 
hoj < 1.0 x 10 3 meV), where fi is the Planck constant di- 
vided by 27r. 

In figure 1, Ss(q,oj)'s at different temperatures arc 
plotted, where we choose the q value corresponding to 
the first-peak position of static structure function (struc- 
ture factor) at each temperature. We recognize a single 
relaxation above 400 K. Here, we call the relaxation as 
the alpha process. At 360 K, a shoulder appears around 
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and the single relaxation at the higher 



temperatures seems to be splitted into two deferent re- 
laxation processes. The relaxation located in the lower 
frequency region is the alpha process. As the tempera- 
ture decreases, the alpha process shifts toward the low to 
region, and the shoulder gradually changes into a peak 
which corresponds to the LEE. The maximum position of 
2.5 x 10 11 s _1 (= ojlee/27t) corresponds roughly to that 
of LEE reported previously in the literature [9j, |l0j. It 
has also been reported that the wlee does not depend 
on temperature generally, while the relaxation time of 
alpha process depends strongly on temperature. Such 
characteristics are well reproduced in Fig 1. 

To know the overall profile of the 5s {q, to), the contour 
maps of log[5s(<z, tL>)] as a function of the logarithm of 
q/lir and to/2w are shown in Fig. 2. At 400 K, a single 
ridge can be recognized from the upper-right corner to 
the lower-left one in Fig. 2(a). As the temperature de- 
creases, the part of the ridge below ui/2tt = 2 x 10 11 s _1 
shifts toward higher q. At the same time, another ridge 
appears almost along the line of lo/2tt = 2 x 10 11 s _1 , 
and hence we can observe that the single ridge at 400 K 
is splitted into two different ones at 300 K. This second 
ridge at to/2n = 2x 10 11 s _1 in Fig. 2(b) corresponds to 
the shoulder observed in Fig.l. At 200 K, the original 
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FIG. 2: Contour maps of log Ss(q, w) as a function of logarithm of to and q at the temperatures of 400 K (a), 300 K (b), 200 
K (c), and 60 K (d). 



ridge below uj/2tt = 2 x 10 11 s" 1 shifts toward higher q 
further, and the second ridge becomes much more dis- 
tinct, so that a peak can be observed at u> = <^lee and 
q/2-K = 5 x 10 9 m — 1 . With decreasing temperature fur- 
ther, the peak position shifts toward higher q, while it 
does not move along the horizontal axis. The peak po- 
sition along the vertical axis corresponds to the spatial 
scale of the LEE at a given temperature. 

In Figure 3(a), the Ss (q, wlee) is plotted as a func- 
tion of the logarithm of g/27r, in order to examine the 
temperature dependence of the spatial scale of the LEE. 
We can notice that all the profiles at the temperatures 
are symmetrical around maxima. The Ss{q, <^lee)'s 
can be well fitted by an empirical function of X(q) = 
A/{(q/q ma x) B + (g/g max )~ B }+constant, where g max , A, 
and B are the position, the height, and the shape factor 
of the peak, respectively. Figure 3 (b) shows the log-log 
plot of the 27r/<7 max vs. the temperature. The 27r/g max 
decreases with decreasing temperature. Non-bonded in- 
ter monomer distance (NBIMD) which is estimated from 
the first peak position of the static structure function also 
decreases from 5.2 x 10~ 10 m to 4.5 x 10~ 10 m with de- 
creasing temperature (not shown). The 2ir/q max at 600 
K is comparable to the the NBIMD. But, it decreases 
further beyond the range of the NBIMD. At 60 K, the 
27r/<7 max is about 1.0 x 10~ 10 m which is smaller than 
a covalent bond length of 1.54 x 10~ 10 m. We notice 
that two straight lines intersect at 230 K corresponding 
to the T g defined in our previous work [24[ , which means 
that the two different power laws hold above and be- 
low T g . After fitting the observed data in Fig. 3(b) to a 
function of 2n/q max — C x T@, we obtained [3 = 0.52 and 
C = 1.1 x 10" u m K-P for 60 K < T < T g , and (3 = 0.97 
and C = 9.3 x 10~ 13 m K _/3 for T g < T < 600 K. 

Here, we discuss the properties of the LEE observed 
in our results. Below T g , the LEE has a major contri- 
bution to the molecular motion, because the character- 
istic time of the alpha process is very large compared to 



that of the LEE. Now, we assume that monomer groups, 
giving rise to the motion associated with the LEE, are 
trapped in potential energy minima. It is expected that 
the monomer groups oscillate around their average posi- 
tions, because large rearrangement of the monomers are 
prohibited below T g . If the displacements of the mass 
centers of the monomer groups from their average posi- 
tions do not correlate and they obey Gaussian distribu- 
tion, the root mean square of the displacement should be 
proportional to T 1 / 2 . Since the the temperature depen- 
dence of the 27r/g max which can be regarded as the aver- 
aged amplitude of the oscillation of the monomer groups 
is almost consistent with the Gaussian approximation, 
we conclude that the motion corresponding to the LEE 
should be described by the above scenario in the temper- 
ature range of 60 K < T < T g , 

As the temperature increases, the characteristic time 
of the alpha process becomes smaller. Then the contri- 
butions of the LEE and the alpha process to the value 
of 27r/g max become superposed in the temperature range 
above T g , so that the temperature dependence of 27r/g max 
becomes strong. We can consider that the monomer 
groups are no longer trapped at fixed positions. The mass 
centers of the groups are drifted by the alpha process. It 
is reasonable that the crossover point of the power law 
appears at T g . 

According to the Gaussian approximation for the dis- 
placements of the mass centers of the monomer groups 
below T g , we estimate the number of the monomers con- 
tained in the group. In this Gaussian approximation, 
square of the averaged amplitude of the oscillation is de- 
rived as (27r/g max ) 2 = 3k B T / m cff w 2 EE , where m cff is 
the effective mass of the monomer groups. Substitut- 
ing 27r/g max = C x T 1 / 2 , the effective mass is written 
as 77i c ff = 3fcB/C 2 a-' 2 EE . From the fitting result of C for 
the data of 60 K < T < T g in Fig 3(b), we obtained 
m c j = 82 g mol -1 which corresponds to 5.9 monomers 
for our model PE. It is smaller than the known character- 
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FIG. 3: Ss(?,wlee) vs. g/27r (a), and 27r/g max vs. temper- 
ature (b). The results of Ss^wlee) from 100 K to 600 K 
are plotted with an interval of 100 K. The maximum position 
of the peak of Ss(q, ojlee) shifts along to the arrow with de- 
creasing temperature. The broken lines in (b) are the fitting 
lines, where the values of the slopes above and below T g is 
0.97 and 0.52 respectively. 

istic sizes concerning to LEE [l8j]. However, our method 
to derive it is much simpler and more direct than the 



others. 

In summary, we have performed the MD simulation 
to investigate the LEE of the model PE system in the 
temperature range above and below T g . We have suc- 
ceeded in observing the LEE at cjlee/27t = 2.5 x 10 11 s _1 
in the contour map of Ss(q,u>). The Ss(9> ^lee) as a 
function of logarithm of q is symmetric around the max- 
imum position of <7 max at a given temperature. The 
27r/g max , which corresponds to the spatial scale of the 
LEE at the low temperatures, depends on temperature 
as 2^/g max ~ T 052 for T < T g and 2^/g max ~ T 97 
for T > T g . Based on the Gaussian approximation for 
the displacements of the monomer groups which give rise 
to the motion regarding the LEE, it is found that the 
number of the monomers included in a group is 5.9. The 
problems of the LEE are not completely solved, but de- 
tailed analysis is under progress to clarify it, which will 
be published in a subsequent paper. 
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